McGill/95-09 
March 1995 



Heavy Resonance Production 
in Ultrarelativistic Nuclear Collisions* 



, - D. SEIBERTt 

■ Department of Physics 

McGill University, Montreal, QC H3A 2T8, CANADA 



in 



> 
in 



X 



ABSTRACT 

Are heavy quarks produced thermally or only by hard parton collisions? What 
is the probability that a produced heavy quark or antiquark is observed in a given 
resonance? 



1. Introduction 



o 

£T) • In this talk, I discuss the problems of heavy quark production and of heavy reso- 

nance production, i.e., are heavy quarks produced thermally or only by hard parton 
, collisions, and what is the probability that a produced heavy quark or antiquark is 

observed in a given resonance? My collaborators are Tanguy Altherr at CERN0, 
who has since died tragically in a climbing accident, and George Fai at Kent State 
University!. 

We use very simple ultrarelativistic nuclear scenarios, basically as proposed long 
ago by Bjorkeni. We assume that the collisions can be divided into the following 
periods: 



1. hard parton collisions, in which essentially all of the final-state entropy is pro- 
duced, for which the typical time scale is about 0.1 fm (we use the standard 
high energy conventions that H — c = ks — 1); 

2. thermal equilibration!, which ends after 0.2 — 0.3 fm, and possibly even sooner!; 

3. expansion of thermally equilibrated gluon plasmai (GP) , or possibly quark gluon 
plasma (QGP); 

4. phase transition, when the temperature T ~ 1 fm -1 ; 

5. resonance gas (RG) expansion; and 

6. freezeout, after which the particles leave the hot matter and flow conveniently 
to the detectors. 

Of course, it is possible that the process of thermal equilibration will coincide with 
some of the later equilibrium processes, but we believe that this will produce only 
small changes to most of our results. It is also possible that there is no phase transition 



but instead some fast but smooth change to the entropy and energy densities as a 
function of T, but again that should not produce large changes in our results. 

We make the following approximations to model the dynamics. First, we assume 
that the hot matter is cylindrically symmetric, so 



f(Px,Py,Pz) = f(PT = \lpl+P 2 y ,Pz)i (1) 

where / is any distribution function. Second, we assume approximate boost-invariance 
as observed in proton collisions^, 

f(r,PT,y)=f(r,p T ). (2) 

Finally, we assume that the hot matter does not expand transversely during the 
collision, so that the volume at proper time r is 

— — = 7tR t, (3) 
dy 

where R is the nuclear radius. 

We also occasionally assume approximate entropy conservation, 

itR 2 ts ~ 3.6 dN/dy, (4) 

where s is the entropy density and dN/dy is the number of produced particles per 
unit rapidity. In that case, we treat the hot matter as an ideal gas with zero chemical 
potential, so 

m _ v 9i f , 3 ^(p) + fc 2 /3^(p) 

Combining eqs. @ and ||, we then estimate T as a function of r for given dN/dy. 

2. Heavy quark production 

In this section I review the theory of heavy quark production. There have been 
three types of calculations of heavy quark production in ultrarelativistic heavy ion 
collisions. The initial production cross sections have been estimated from perturbative 
quantum chromodynamics (QCD), either by simply taking cross sections for proton 
collisions! and scaling by A 2 , where A is the number of nucleons in the nucleus, or 
in a more sophisticated manner by using the nuclear (instead of the bare nucleon) 
parton structure functionsEUl. The results obtained with these two methods do not 
differ very much. 

Thermal production cross sections have been calculated by taking the production 
matrix elements, convoluting with the four-volume element and the thermal distri- 
bution functions, and integrating over assumed collision historiesifflifl. These 
thermal calculations have grown in sophistication with time, but the results have not 
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changed very much. Our major addition to these calculations has been the inclusion 
of heavy quark production through thermal gluon decay, which has been previously 
neglected but is the dominant term in the weak coupling due to the anoma- 

lously large thermal gluon widthEl. For typical values of the strong coupling constant 
probed in nuclear collisions, this new production term is comparable to those previ- 
ously calculated, although it does not dominateB. 

Finally, there is the parton cascade modeli@'0, which attempts to include both of 
the previous calculations by using a perturbative QCD cascade, in which interactions 
are more or less arbitrarily cut off at some lower momentum so that cross sections 
remain finite. This model seems to work reasonably well for strange (s) quark pro- 
duction, which at least seems to be correctly predicted for proton collisions and scales 
with projectile and target in a reasonable manner. Results for charm (c) and bottom 
(b) production are probably not reliable, since proton collision results are not repro- 
duced and the cross sections scale as A 5 / 3 in proton-nucleus collisions, instead of A 
as observed by experimenters. 

In my opinion, the best density estimates and dominant quark production mech- 
anisms for an Au+Au collision at RHIC energy (^/s = 200 GeV/nucleon) are: 

s: Pure thermal production gives dN s /dy ~ 100, the parton cascade model gives 

dN s /dy ~ 50, so I expect that dN s /dy ~ 100, mostly from thermal collisions, 
c: Thermal production gives dN c /dy ~ 1, perturbative QCD gives dN c /dy ~ 1, 

so I expect that dN c /dy ~ 1 with roughly equal production from hard parton 

collisions and thermal collisions. 
b: Thermal production gives dN b /dy < 10~ 3 , perturbative QCD gives dN b /dy ~ 

0.02, so I expect dN^/dy ~ 0.02 mostly from hard parton collisions. 



3. Freezeout conditions 

Once we have our heavy quarks, the next question is, "When do they freeze out?" 
We attempted to solve this problem with a simple model of heavy quark production 
and dynamics in an ultrarelativistic nuclear collisionll. We first assume that all heavy 
quark- antiquark pairs are produced at time r = 0; this is reasonable for c and b quarks 
even if they are produced thermally, since thermal production is also concentrated 
at early times, although it is not such a good assumption for s quarks. We then 
assume that these quarks quickly thermalize, and that their subsequent trajectory is 
a random walk in a thermal bath, with collision frequencies taken from Pi S ar s ki§; 
the bath temperature is estimated by assuming entropy conservation. Finally, we say 
that the quarks have frozen out when either (i) their mean free path is larger than 
the mean distance to leave the hot matter, or (ii) their random walk has taken them 
out of the hot matter. 

This is of course complicated by the fact that we expect a phase transition to RG 
during the evolution of the hot matter. We take this into account by assuming that 
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A dN/dy Q T f /T c r f (fm/c) r Q (fm) n f 
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Table 1: Frcczeout parameters, adapted from Ref. 2. 

the interaction cross section per unit entropy is approximately the same in the two 
phases, and that in the RG phase the quarks are contained in mesons which we model 
as qq pairs to estimate the mesonic mean free paths. We thus obtain different results 
by varying the transition temperature, T c , and the ratio of the number of degrees of 
freedom in the two phases, v. 

I show results in Table 1 for S+S collisions at SPS energy (-y/s = 20 GeV/nucleon), 
and for Au+Au collisions at SPS, RHIC, and LHC (^/s = 7 TeV/nucleon) energies, 
averaged over T c = 150 and 200 MeV, and over v — 5 and 10. Here Tf is the freezeout 
temperature, tj is the proper time, rq is the distance the quark has moved from its 
point of creation, and n/ is the mean number of collisions before freezeout. This last 
is the most important number - when rif ~ 1 (S+S collisions at SPS), statistical 
recombination models should not be expected to work very well, while when nf ^> 1 
(Au+Au collisions at all energies considered) we expect that statistical models should 
describe heavy resonance data very well. 

4. Statistical recombination 

Our statistical recombination model! is very similar to those used by recent 
author JllllIB. The main difference is that, in addition to ensuring that the quark 
and antiquark densities, respectively pq and pg, are conserved at freezeout, we en- 
sure that the local quark-antiquark density, p^(x, x), is also conserved. This is not 

(2) 

trivial, as typically p^(x, x) ^> PqPq- We thus include three chemical potentials: 
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1. quark, 

2. antiquark, fig ^ — fig, since the heavy quarks are not in chemical equilibrium; 

3. pair, Hqq- 

The chemical potential for resonance i is 

tk = E ( k P»Q + k ?^Q + • ( 6 ) 

Here fc t - ^ (&f^) is the number of quarks (antiquarks) of flavor Q, and k\®^ is the 

number of pairs (the smaller of kf^ and fcj ). We include all confirmed meson and 
baryon resonances^ in our statistical recombination model. 

We estimate the two-particle density (and thus Hqq) from the freezeout conditions 
obtained from our simulation. 

Results are shown in Table 2, again averaged over T c and v. Although the corrections 
to s resonance production from the pair chemical potential are less than order unity, 
the corrections to c and b resonance production are huge, giving orders of magnitude 
changes in fugacities. Thus, this pair chemical potential (or its equivalent) must be 
included to calculate c and b resonance production. 

Finally, this recombination model can be used to gain insight into predicted phe- 
nomena such as suppression of J/if> resonance production in ultrarelativistic nuclear 
collisions^!. Predictions for the fractions of cc and bb pairs that freeze out as cc and bb 
mesons are given in Tables 3 and 4. These quantities serve as surrogates for the J/ip 
and T suppression, as these resonances are the most common sources for observed 
J/0 and T resonances. It is worth noting that the heavy mesons are all predicted 
to freeze out in the RG in ultrarelativistic nuclear collisions, so their suppression 
probably depends most strongly on (and thus carries the most information about) 
their interactions with the RG, and not with the QGP; thus, their dynamics is only 
weakly dependent on the presence or absence of QGP. Otherwise, it is obvious from 
the strong dependences on T c and v that the theoretical uncertainty in these calcu- 
lations is huge, so much work remains to be done on the problem of how produced 
heavy quarks emerge as observed heavy resonances. 



5. Conclusions 

We have calculated heavy quark production; s quarks are mostly produced ther- 
mally, b quarks are mostly produced by hard parton collisions, whil c quarks are 
produced by both mechanisms. We have attempted to estimate when freezeout will 
occur for the various heavy quark species, and whether there are enough collisions 
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Tabic 2: Freezeout densities, adapted from Ref. 2. 



T c (MeV): 150 200 

A dN/dy v. 5 10 5 10 

132 85 0.066 0.075 0.029 0.021 

197 1000 0.021 0.002 0.008 0.001 

197 2000 0.009 0.004 0.004 0.001 

197 3500 0.009 0.010 0.005 0.002 



Table 3: cc meson fractions, from Ref. 2. 
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Table 4: bb meson fractions, from Ref. 2. 
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that statistical recombination models should be reliable guides to heavy resonance 
production; the latter seems to be true for Au+Au collisions at SPS energy and 
above, but not necessarily for S+S collisions at SPS. We found that heavy quarks 
typically freeze out from the RG, and not from the QGP, so the relative numbers of 
the various heavy resonances are only weakly dependent on the presence or absence 
of QGP. Finally, we found that it is necessary to include a pair chemical potential in 
statistical recombination models, or to otherwise allow for the fact that heavy quarks 
and antiquarks are always produced in coincidence. 

This work was supported in part by the Natural Sciences and Engineering Research 
Council of Canada, and in part by the FCAR fund of the Quebec government. It 
is based on earlier work supported in part by the U.S. Department of Energy under 
Grant No. DOE/DE-FG02-86ER-40251, and in part by the North Atlantic Treaty 
Organization under a Grant awarded in 1991. 
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